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1. Introduction 

The Astrophysical Multimessenger Observatory Network (AMON) links existing and forth¬ 
coming high-energy astrophysical observatories into a single virtual system, capable of sifting 
through the various data streams in near real-time, identifying candidate and high-significance 
multimessenger transient events, and providing alerts to interested observers [1, 2]. The AMON 
framework also enables searches for coincidences in archival data. This helps us to better compre¬ 
hend the datasets and also to explore different statistical approaches to generate realtime AMON 
alerts for the network’s followup partners. 

Many theoretical models predict that detectable neutrino emission will be paired with a prompt 
electromagnetic signal [3,4]. Recently, the IceCube Observatory reported an excess of astrophysi¬ 
cal neutrinos in the eV energy scale [5]. These motivate us to perform archival coincidence 

analyses between the IceCube public data and the Fermi LAT public photon datasets. Our goal is 
to look for statistically significant coincidences between high-energy neutrinos and gamma-ray 
signals to find jointly-emitting neutrino/gamma-ray transient sources. 

This paper is organized as follows: The detail of the datasets we use is discussed in section 2. 
Different statistical methods are discussed in section 3. Section 4 shows the results of the analyses. 
Conclusions and future work are presented in section 5. 

2. Datasets 

This analysis has been performed using IceCube and Fermi LAT public data on a period of 
temporal overlap between the two observatories. IceCube public data includes the 40-string and 
59-string configurations of IceCube observatory (IC40 and IC59). IC40 dataset starts on April 6, 
2008 and ends on May 20, 2009 and IC59 dataset is between May 20, 2009 and May 31, 2010. The 
public Fermi LAT dataset starts on August 4, 2008 and has been continuously available since then. 
The overlap period of IC40 and Fermi LAT data is therefore between weeks 9 and 50 of the Fermi 
LAT data. This period for IC59 is between weeks 50 and 104. 

The IC40 public data only contains neutrinos from the northern hemisphere (positive declina¬ 
tion). In total, there are about 14,000 neutrinos in this dataset. For the purpose of this analysis, we 
only consider the northern hemisphere data of IC59 as well which contains about 43,000 neutrinos. 
The Fermi LAT data that is being used is selected by imposing additional criteria using Fermi LAT 
purest (i.e. lowest instrumental background) analysis class (Pass 7-V6 Ultraclean). In this class, 
the photon events that are detected while the telescope is repositioning, or is in pointing mode, 
or pointing to the Sun are removed. In addition, only photon events with spacecraft zenith angle 
smaller than 65° and energies above 200 MeV are accepted. 

Applying all these criteria, the dataset reduces to approximately 4.1 million photon events 
for the period overlapping IC40 and 5.5 million photon events for the period overlapping IC59. 
Figure 1 shows the Fermi LAT exposure corrected map for the IC40-I-IC59 period (on the left) 
and the observed neutrinos in northern hemisphere of IC40-I-IC59 data (on the right). There is an 
excess flux above the horizon in the neutrino map that is presumably the result of cosmic ray muons 
that are incorrectly reconstructed as coming from the northern sky. Both maps are in equatorial 
coordinates. 
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(a) Fermi LAT (b) IC40+IC59 


Figure 1: Sky map in equatorial coordinates for (a) Fermi LAT photons on a period of temporal overlap 
with IC40+IC59 and (b) IC40+IC59 neutrinos. 


To start the analysis, we calculate the angular separation between each photon event and each 
neutrino and identify all pairs with smaller than 10° separation and |At| < 50 s. This value of 
angular separation is chosen due to the large angular uncertainty of LAT data at lower energy. 
IceCube neutrino angular uncertainty is about 1°. A 10° angular separation between a photon- 
neutrino pair is an upper limit that guarantees the doublet searches. We also apply a temporal cut 
(AT) of ±50 s between neutrino and photon arrival times. 

3. Statistical Methods 

3.1 Likelihood 

To perform the archival coincidence analysis between IC40/59 and Fermi LAT data, an un¬ 
binned log-likelihood function (A) is used. This metric depends on event arrival direction and its 
uncertainty as well as a background rejection term and is defined as follows: 

A = 21n (/iAT(x|xy)P[c(x|xv)) - 2In {B{xj)) , (3.1) 

where A is the likelihood parameter, Xy and Xy are the arrival direction of the photon and neutrino 
events, respectively, x is the best fit position, B{xy) is the Fermi LAT background rejection term that 
will be briefly discussed in Section 3.2, Flat(^|-^ 7 ) and Pic{x\xy) are the energy dependent point 
spread functions (PSF) of the LAT and IceCube observatories, respectively. In this metric, large 
values of A indicate a higher probability correlated pair. 

The PSF for IceCube is given by a Gaussian function with an energy dependent spread [6]. 
The PSF for Fermi LAT is given by a two parameter King function [7], deriving its energy depen¬ 
dency from empirically determined values. The best fit position of any doublet, x, can be obtained 
analytically by solving a cubic equation. 

3.2 Background Rejection 

The Fermi LAT background rejection term is a distribution that is proportional to the event 
rate over the exposure. For photons that end up in the galactic plane, the value of B{x) is large, 
and therefore the pair will be given a less significant weight. A HEALPix map [8] with nside=64 
(mean spacing of 0.91°) is used to construct the background rejection map. 
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3.3 Null and Signal Distributions 

We generate a series of 10,000 scrambled datasets to estimate the null distribution. We scram¬ 
ble time and right ascension of each neutrino event while keeping declination and intrinsic neutrino 
properties (such as energy and directional resolution) unchanged. Time is shuffled using a uniform 
random number from a total time window of overlapping neutrino and photon datasets. We note 
that since LAT’s motion is complex, we do not time-scramble the photon events. Applying the 
spatial and temporal cuts on all of the IC40 and Fermi LAT neutrino-photon (v — 7 ) pairs using the 
method discussed in Section 2, on average 2207±40 pairs in each of the 10,000 null datasets are 
left. Then a log-likelihood (A) value is calculated for each v — / pair from the Equation 3.1. 

Additionally, a series of 10,000 signal tests are created by injecting a limited number of forced 
coincidences into the null data. Each signal dataset is obtained by matching a LAT photon event 
with each of the IceCube neutrino events. First, a random uniformly distributed source position is 
chosen on the sky as the true event position around which the messenger events will be distributed. 
Then, the neutrino and photon events are randomly given a position on the sky centered around the 
north pole being the true event position. This way, 0 is uniformly distributed and 6 is drawn from 
the PSF of each messenger. These random positions are then rotated into space by transforming 
the north pole into the randomly chosen true event position. The photon energy is drawn from the 
convolution between the normalized power law distribution with an arbitrary spectral index of 2.2 
and the LAT exposure function. Temporal coincidence is ensured through small time stamp offset. 
The angular separation between each photon and each neutrino is calculated and only pairs within 
the cut of 10 ° are accepted and then A is calculated for each pair. 

Figure 2 shows the normalized cumulative histograms of A values for both null and signal 
distributions using the described 10,000 datasets for IC40 and Fermi LAT coincidences. 



Figure 2: Histograms of log-likelihood (A) values of the null and signal distributions. Both distributions are 
normalized and cumulative. Signal distribution shows higher A values, as expected. 
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3.4 Defining Statistical Excess 

We perform an Anderson-Darling (AD) test [9] on the null and signal likelihood distributions 
to determine if there is a statistical excess. A p-value can be obtained from this test statistic using 
an AD k-sample test [10]. We inject signal photons randomly chosen from a single signal run with 
spectral index of 2.2. The results of obtained AD test statistic p-value is plotted versus the number 
of injected signal photons (A^sig) in Figure 3. The injected signals are added to the null datasets 
and the new signal datasets is tested against the null hypothesis. This plot shows the ability of this 
method to distinguish injected signal from background. In other words, it helps us to obtain an 
estimate for the analysis sensitivity. We conclude from this plot that our analysis is sensitive to 
presence of < 100 v — 7 pairs in full dataset. 
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Figure 3: Anderson Darling (AD) test statistic p-value versus number of injected signal photons (A^sig)- This 
plot helps us to obtain an estimate for the analysis effectiveness. The p-value obtained from the AD test 
suggests about 70-100 signal events. 


4. Results 

4.1 Unblinding Data 

Here we only present the results of unblinding the IC40 and Fermi LAT data. Unblinding 
the IC40 data and performing the coincidence analysis with Fermi LAT data result in 2138 co¬ 
incidences applying the cuts of A0 < 10° and |At| < 50 s. The AD test statistic on this dataset 
provides a p-value of about 4%. This p-value is consistent with simulations including the injection 
of 70-100 signal events. The A distribution plot for unblinded data is presented in Figure 4. The 
null hypothesis and a theoretical signal model with 70 injected signals {N^ig=10) are also plotted for 
comparison. Note that all distributions are normalized and multiplied by 2138, which is the total 
number of coincidences after the unblinding. The bottom plot in Figure 4 shows the residuals of 
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data and A^sig=70 versus the null hypothesis to better represent the differences between the distribu¬ 
tions. Adding 70 injected V — 7 pairs to the null data increases the number of high-A events. This 
explains the decrease observed in the residual plot of the two normalized cumulative distributions 
(A^sig=70 and the null hypothesis) shown in red dashed line. Although the cumulative distribution 
is not a particularly good match to the simulated distribution for A^sig=70, we nonetheless take this 
as our model in defining a future three vetting tests for the data. 
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Figure 4: Top: cumulative distribution of log-likelihood values of the IC40 and Fermi LAT coincidences 
after unblinding the data. A theoretical signal model with 70 injected signals (A^sig=70) and the null hypoth¬ 
esis are plotted for comparison. All distributions are normalized and multiplied by 2138, which is the total 
number of coincidences after unblinding the data. Bottom: The residuals of data and Asig=70 versus the null 
hypothesis. 


4.2 Vetting the Signal 

Three different further tests are performed to explore whether real cosmic V — 7 pairs are 
present in the data at Vsig=70 level. We perform the tests only on events with high-A values which 
are more likely to be signal. To find the threshold above where we consider the V — 7 pairs (Acut), 
we calculate the signal to noise ratio (SNR) for different Acut values. SNR is defined by the number 
of injected signals with A > Acut divided by the square root of total number of pairs (injected signal 
+ null data) also with A > Acut- We find the maximum SNR at Acut=l 1- Only events with A > Acut 
are considered in the following tests. 

The first test is multiplicity, i.e. the number of photons in coincidence with a single neutrino. 
The mean multiplicity of events with high-A is 2.17 for real data whereas this number is 2.08 
for the null hypothesis with a standard deviation of 0.15 between the corresponding values of the 
10,000 datasets. Comparison between the mean multiplicity values of V — 7 pair coincidences in 
data and the null hypothesis indicates no significant signal excess above the background expecta¬ 
tions with a p-value of 0.25. 
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(a) At (b) Clustering 

Figure 5: Vetting the signal: (a) At: the histogram of the time difference between the neutrino and the 
photon in each pair is plotted and compared with a flat uniform distribution, (b) Clustering: source map of 
the high log-likelihood v — / pairs (shown by grey circles). In total, there are six v — / pairs that lie within 
2° of another pair, which is less than the result we get from null hypothesis. Red crosses are the locations 
of the four neutrinos that contribute to the total of six v — / pairs. These tests indicate no signiflcant signal 
excess in IC40 and Fermi LAT datasets. 


The next test is to look for a significant difference between different At bins. At is the time 
difference between the IC40 neutrino in coincidence with the Fermi LAT photon event. Figure 5(a) 
shows the At histogram for five time bins compared to a uniform distribution. The reduced chi- 
square value is approximately 0.5, which indicates consistency with a flat uniform distribution, 
further suggesting absence of signal. 

Finally, we test clustering of the high log-likelihood V — / pairs which would reveal bright/repeating 
sources. Figure 5(b) shows the source map of these pairs. We see that only six v — / pairs lie within 
2° of another pair. In total, four neutrinos are contributed in creating these six pairs which are in¬ 
dicated by red crosses on the sky plot. Six v — / pairs are shown to be less than the average result 
of 12.9 such clustered V — / pairs we get from our 10,000 scrambled datasets; roughly no evidence 
for a contribution for cosmic V — / pairs. 

5. Conclusions 

We performed an archival analysis on neutrinos from IceCube observatory in coincidence with 
Fermi LAT photon events, both from public datasets. Several statistical tests on observed data using 
the background and signal datasets were conducted. The Anderson-Darling test statistic showed 
about 70 signal out of 2138 found coincidences in IC40-Fermi LAT analysis, however multiplicity, 

AT, and clustering tests showed no significant signal excess. The results of IC59 and Fermi LAT 
and the combined IC40/59 datasets and Fermi LAT will be presented in a future publication. 
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